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I.  INTRODUCTION 


A  problem  of  continued  interest  is  that  of  antennas  mounted  on  finite 
conducting  bodies.  Antennas  on  spherical  bodies  have  been  analyzed  [1], 
as  well  as  antennas  on  more  general  structures [2,3] .  The  authors  have 
previously[4]  presented  a  surface  patch  moment  method  solution  for  geom¬ 
etries  involving  wires  and  rectangular  plates,  including  wire/plate  and 
plate/plate  junctions.  In  that  work,  the  wire  attachment  point  was  re¬ 
stricted  to  be  O.lX  or  more  from  any  plate  edge.  This  report  presents 
work  which  allows  attachment  points  to  be  near  or  at  a  half-plane-type 
edge,  and  also  serves  as  a  start  to  the  problem  of  wire  attachments  near 
the  edge  of  two  plates  intersecting  at  an  arbitrary  angle  (a  wedge). 

Also  presented  are  canonical  solutions  to  the  problem  of  a  monopole  near 
or  at  the  edge  of  a  perfectly  conducting  half-plane.  These  solutions  are 
in  terms  of  near  field  parameters,  i.e.,  currents  and  input  impedance. 
Far-field  results  for  antennas  near  half-plane  edges  are  available  in  [51. 

Various  types  of  solutions  are  presented  for  each  of  two  related  geom¬ 
etries:  a  vertical  monopole  near  an  edge  and  a  radial  monopole  mounted 
on  an  edge.  The  first  type  of  solution  is  a  rigorous  hybrid  Green's  func¬ 
tion  moment  method  solution  for  an  antenna  near  or  at  the  edge  of  a  per¬ 
fectly  conducting  half-plane.  This  solution  uses  an  integral  expression 
for  the  exact  half-plane  Green's  function,  with  the  source  and  image  terms 
separated  from  the  edge  diffracted  field.  This  procedure  serves  as  an 
accurate  canonical  solution  for  comparison  with  other  solutions.  Other 
results  of  this  solution  are  simple  closed  form  expressions  for  the  near 
zone  surface  currents  for  a  source  near  or  at  an  edge,  which  are  later 
used  to  synthesize  a  surface  patch  moment  method  attachment  mode.  Dif¬ 
ferent  Green's  functions  are  used  depending  on  whether  the  antenna  is  near 
the  edge  or  on  the  edge. 

The  second  type  of  solution  presented  is  similar  to  the  first, 
except  the  GTD  is  used  to  find  the  diffracted  field,  rather  than  using 
the  exact  Green's  function.  This  solution  is  essentially  the  MM/GTD  hybrid 
method  developed  by  Thiele  and  Newhouse[2].  Agreement  with  the  MM/Green's 
function  is  good  for  vertical  monopoles  as  close  as  .2X  from  the  edge. 


Since  GTD  is  not  accurate  for  source  and  field  points  both  very  near  the 
edge,  the  MM/6TD  solution  was  not  carried  out  for  radial  monopoles  mounted 
on  an  edge. 

« 

The  third  solution  involved  constructing  a  new  expansion  mode  for 
wires  attached  near  or  at  an  edge,  for  use  in  the  surface  patch  MM  solu¬ 
tion  [4],  The  new  mode  was  based  on  the  form  of  the  surface  currents  found 
via  the  exact  Green's  function.  Also,  a  special  mode  to  enforce  the  proper 
edge  singularity  for  the  component  of  surface  current  parallel  to  the  edge 
was  developed.  This  mode  was  found  to  have  practically  no  effect  on  the 
input  impedance. 

II.  SOLUTION  FOR  ANTENNA  NEAR  EDGE  OF  HALF-PLANE 

USING  EXACT  GREEN'S  FUNCTION 

The  problem  of  a  wire  antenna  in  free  space  can  be  solved  via  the 
moment  method  solution  of  an  integral  equation  using  the  free  space  Green's 
functionlB).  If  the  wire  antenna  is  in  the  presence  of  another  scattering 
structure  there  are  two  possible  modifications.  The  first  is  to  replace 
the  free-space  Green's  function  by  the  special  Green's  function  for  the 
particular  scattering  structure.  The  second  approach  is  to  retain  the 
free-space  Green's  function,  but  suitably  expand  the  currents  on  the  wire 
and  the  scattering  structure,  and  solve  for  the  expansion  coefficients 
in  the  MM  solution.  The  advantage  of  the  first  approach  is  that  no  new 
unknowns  are  introduced  into  the  MM  solution  and  the  size  of  the  impedance 
matrix  does  not  increase  over  that  of  the  isolated  wire,  while  the  dis¬ 
advantage  is  that  the  Green's  function  of  the  scatterer  must  be  known  in 
a  computationally  efficient  form.  The  advantage  of  the  second  approach 
is  that  the  Green's  function  of  the  scatterer  need  not  be  known,  while 
the  disadvantage  is  that  new  unknowns  are  introduced  and  the  size  of  the 
required  impedance  matrix  can  increase  substantially  if  the  scatterer  is 
electrically  large. 

The  first  approach  is  used  in  this  section  and  in  Sections  III  and 
IV,  while  the  second  approach  is  used  in  Section  V. 
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Figure  1.  Geometry  and  coordinate  systems  for  a  vertical  monopole 
near  the  edge  of  a  half-plane. 


A.  Method  of  Solution 


The  geometry  of  the  problem  is  shown  in  Figure  1.  The  wire  and  half¬ 
plane  are  for  simplicity  considered  to  be  perfectly  conducting.  The  basic 
idea  of  the  solution  is  to  expand  the  wire  currents  using  piecewise-sinus- 
oidal  modes,  find  the  fields  of  test  modes  (also  PWS)  in  the  presence  of 
the  half-plane,  and  compute  an  impedance  matrix: 

^i 


3 


I 


where  current, 

Ii(y)  is  the  j-th  expansion  mode  current, 
and  Eyy(y|yQ)  is  the  half-plane  Green's  function  representing  the 

Ey  field  due  to  a  y  directed  infinitesimal  electric  dipole  of 
unit  moment. 

Since  IJy)  =  I -(y)  we  have  a  Galerkin  solution,  and  Z. .  =  Z-^.  The  wire 
currents  and  input  impedance  then  follow  from  the  solution  of  a  matrix 
equation  of  the  form  [ZlII]  =  [Vl,  where  (Zl  is  the  impedance  matrix,  [1] 
is  the  column  vector  of  current  expansion  coefficients,  and  [V]  is  the 
vector  of  mode  voltages. 

Thus,  this  solution  depends  upon  the  availability  of  the  half-plane 
Green's  function  in  a  form  suitable  for  fast  and  accurate  evaluation.  It 
was  first  attempted  to  use  the  familiar  form  of  the  half-plane  Green's, 
function  expressed  as  an  eigenfunction  expansion,  such  as  given  by  Tai 
f5).  This  attempt  failed  for  various  reasons. 

First,  the  eigenfunction  expansion  of  the  half-plane  Green's  function 
consists  of  an  infinite  series  of  infinite  integrals,  making  evaluation 
in  the  near  zone  very  difficult  due  to  slow  convergence.  The  slow  conver¬ 
gence  is  due  to  the  fact  that  the  eigenfunction  expansion  contains  the 
singularities  for  the  source  and  image  terms.  It  was  found  that  the  source 
(or  incident)  terms  could  be  subtracted  off,  and  expressed  in  simple  closed 
form,  however,  the  authors  were  not  able  to  do  the  same  for  the  image  (or 
reflected)  field.  Thus,  the  evaluation  of  self  or  overlapping  mutual  im¬ 
pedance  terms  would  have  been  very  difficult  to  do  accurately.  Also,  the 
eigenfunction  solution  did  not  lend  itself  to  small  argument  approximations 
for  evaluation  of  near  zone  surface  currents  in  simple  closed  form.  It 
was  found  that  the  first  term  of  the  infinite  series  would  give  the  correct 
edge  behavior,  but  not  the  correct  source  singularity;  the  entire  series 
was  needed  for  that.  For  these  reasons  another  type  of  solution  for  the 
half-plane  Green's  function  was  sought. 
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In  1953,  Senior  f71  gave  the  solution  for  diffraction  of  a  vertical 

electric  dipole  by  a  half-plane.  His  solution  was  expressed  in  terms  of 

derivatives  of  the  scalar  diffraction  solution,  which  was  expressed  in 

integral  form.  As  will  be  seen,  this  solution  overcomes  the  problems  which 

were  met  with  the  eigenfunction  solution. 

I 

B.  Formulation  of  MM/Green's  Function  Solution 


For  a  vertical  infinitesimal  electric  dipole  of  moment  Ij,  in  the  pres¬ 
ence  of  a  conducting  half-plane,  the  total  field  can  be  written  as  (71 
(for  e"^^  time  dependence): 


y  4nj 


cos  ^  cos  -x— 


(2) 
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The  other  field  components  (E  ,E  ,H  )  are  not  needed  in  this  study 

X  z  y 

but  can  be  obtained  from  [7l. 


The  incident  and  reflected  electric  fields  can  be  written  in  terms 


of  a  Hertz  potential: 


111 

-me  _  _  _o 
■y  “  4ifj 
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e-JkR 

TIT" 


y  4itj  V9y 
Usin^  the  relation  [7], 


;  Hj^^kR  cosh  u)du  =  2j  , 


in  (5)  gives 


-jkR  <*>  iy\ 

Ir  =  2j  ^ - ;  Hj'^^kR  cosh  u)du, 
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■gjT  3^  ♦  written  as 


^  1  ^  ^ ^  *(i7  *  -(^0  ■ '') 
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where  we  have  defined 


=  ^  /  Hj^^(kR  cosh  u)du; 
^‘R 


(10a) 


=  ^  /  Hp^(kR'  cosh  u)du 

Ur, 


(10b) 


Note  that,  using  (7),  and  can  be  expressed  in  terms  of  a  finite 
integral : 


dR  1  i  (2) 

^ ^  H|^^(kR  cosh  u)du. 


(11a) 


dR'  1  j  1^'  (2), 

"  "  7  — RT*”  “  2  ^  Hj^^(kR'  cosh  u)du. 

It  can  be  shown  that  and  are  not  singular  at  R=0  or  R'=0. 


(11b) 


The  geometrical  optics  (incident  and  reflected)  fields  can  now  be 
easily  recognized,  via  (6),  and  separated  from  the  diffracted  field: 


rinc  ^  nref 

y  y 


-  rg*o.  ^  cdif 

■  y  y 


(12) 


E  = 
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+  E 


dif 
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The  diffracted  electric  field  is  now  in  a  form  which  has  no  singularities 
at  the  source  or  image. 


Since  the  geometrical  optics  terms  have  been  separated  from  the  dif¬ 
fracted  field,  the  impedance  matrix  [Z]  can,  for  computational  purposes, 
be  broken  into  two  parts;  [Z]  =  [Z^]  +  [AZ],  where  [Z^] represents  the  im¬ 
pedance  matrix  of  a  monopole  on  an  infinite  ground  plane  contri- 

d  i  f  ^ 

bution)  and  (AZ  ]  represents  the  E^  field  contribution  which  accounts 

for  the  edge.  If  PWS  modes  are  used  for  expansion  and  testing  the  [Z^l 

elements  can  be  computed  in  closed  form  using  Richmond's  subroutines  [8], 
eliminating  the  difficulties  which  would  be  encountered  if  this  evaluation 
were  done  numerically.  This  is  an  important  and  useful  distinction  from 
other  similar  solutionsll] .  The  remaining  task  is  to  evaluate  [AZ] . 


The  PWS  modes  are  arranged  on  the  wire  as  shown  in  Figure  2a,  num¬ 
bered  from  top  to  bottom.  The  test  modes  are  positioned  on  the  wire  axis 
(at  x^,  Zg=0)  while  the  expansion  modes  are  placed  on  the  wire  surface 
at  (x  ,z  =a),  where  a  is  the  wire  radius.  Note  that  the  modes  are  PWS 
dipoles  except  for  the  bottom  mode,  which  is  a  PWS  monopole  with  its  ter¬ 
minal  at  the  base  of  the  wire.  This  allows  non-zero  current  at  the  attach¬ 
ment  point.  Figure  2b  shows  a  typical  test  and  expansion  monopole,  with 
coordinates,  which  make  up  the  PWS  modes.  All  segments  are  chosen  with 
the  same  length,  d  =  =  It^-tj].  The  current  on  the  monopoles  shown 

in  Figure  2b  can  be  written  as 
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Figure  2a.  PWS  expansion  mode  arrangement  on  wire 


Figure  2b.  Typical  test  and  expansion  monopoles  and  coordinates. 


for  ^2^1* 


sin  k(t2-yQ) 
sin  kd 


sin  k(ep-y) 

'e'y’  ■  — Wkd  '<”•  «2=^1'  *<-‘>  V=l- 

The  contribution  to  the  mutual  impedance  between  these  two  monopoles  due 
to  the  diffracted  field  can  then  be  written  as 

t  ®2 

AZt.e  =  -  ^(^0^  J"  ^e^y^  Ey'^^(y|yo)t^y  % 

where  (y|yo^  given  by  (13).  Fortunately,  the  double  differentiation 

in  (13)  can  be  removed  by  integration  by  parts,  giving 


4  ®2  m. 


AZ 


sin  k(e2-y)sin  k(t2-yQ) 


2 

sin  kd 


=  -  f  r  0 

t,e  ^  4Tij 

yo=ti  y=ej 

-ik  (b  ***0  f2l 

.  cos  I  cos  -j  Hj'^^(kRj) - 2' •  •  "  -  ^ 


sin  kd 


k  cos  k(e2-t2-y+yQ) 

- —x - n 

sin  kd 


dy  dy^ 


(15) 


The  expression  (15)  has  been  simplified  considerably  due  to  the  fact  that 
dR  dR  * 

Ti  (yUg)  =  ^  ^y|yo^  ^o"®’  because  PWS  modes  with  a  free- 

space  wavenumber  were  used.  Also,  (15)  does  not  contain  certain  terms 
such  as  end  charges  which  would  cancel  when  two  monopoles  are  combined 
to  form  a  dipole.  Thus,  [AZl  ran  be  calculated  with  two  numerical  inte¬ 
grations  for  (15)  and  one  for  (11).  These  integrations  are  well  behaved, 
typically  requiring  only  10-20  Simpson  intervals  each.  This  ease  of  com¬ 
putation  of  [AZ]  is  a  direct  result  of  the  fact  that  the  incident  and  re¬ 
flected  field  components  were  separated  from  the  total  field,  allowing 


lu 


■J 


the  nearly  singular  portion,  [Z^jj  ,  to  be  evaluated  in  closed  form.  Also, 
Equation  (15)  is  valid  for  monopole  orientations  where  ej^>e2  or  tj>t2, 
which  are  required  for  the  complete  dipole  impedances. 

C.  Results 

• - 

The  solution  described  above  gave  a  quickly  converging  result  for 
the  wire  currents,  as  shown  in  Table  I  below.  This  table  gives  AZ  =  Z^-^- 
Zq,  where  Z^^  is  the  input  impedance  of  a  X/4,  a= .OOlX  monopole  near  the 
edge  of  a  half-plane,  and  Z^  is  the  input  impedance  of  the  same  antenna 
mounted  on  an  infinite  ground  plane,  is  the  distance  of  the  antenna 
from  the  edge.  Results  are  shown  for  1,  2,  and  3  modes  on  the  monopole. 

A  delta  gap  generator  at  the  base  of  the  monopole  was  used. 

Table  I 


Xo/X 


AZ,  1  mode 


AZ,  2  modes 


AZ,  3  modes 


From  the  above  data  it  is  seen  that  the  solution  has  practically  con¬ 
verged  with  two  wire  modes.  AZ  vs  is  plotted  in  Figure  7  in  Section 
1 1  IB,  where  it  is  compared  with  the  solution  discussed  in  that  section. 

In  Section  VB  vs  is  plotted  in  Figures  11  and  12  and  compared  with 
the  results  of  that  section.  Computation  time  for  this  solution  iS| about 
30  seconds*  when  2  expansion  modes  are  used. 

For  Xq  =  .05  with  two  wire  modes  the  impedance  matrix  elements  are 
(23.89-j410.5)  (12.64+j3l8.5) 

I^ol  = 

.(12.64+j318.5)  (6.71-j224.1) 

(2.51-j4.12)  (4.25-j2.38)  ' 

(4.25-j2.38)  (4.85-jl.67)  J  . 

The  wire  modes  are  numbered  from  top  to  bottom  (attachment  point). 

The  solution  presented  in  this  section  is  valid  for  any  distance, 

Xq,  from  the  half-plane  edge  except  when  the  antenna  is  within  a  fel#  wire 
radii  of  the  edge.  In  this  case  numerical  problems  are  encountered  since 
the  fields  contain  the  \/Jx  edge  singularity.  The  solution  presented  in 
Section  IV  overcomes  this  difficulty.  *■ 

D.  ({>-variation  of  Wire  Current 

The  solution  in  the  previous  section  used  an  axial  current  filament 
for  testing  and  a  circumferentially  uniform  surface  current  located  on 
the  surface  of  the  wire  for  expansion.  Thus,  the  solution  did  not  account 
for  or  give  any  information  about  the  circumferential  variation  of  the 

*A11  times  reported  here  are  for  Datacraft  6024  -  about  8  times  slower 
than  IBM  360/165. 


current  around  the  wire  surface.  This  aspect  of  the  solution  was  of  in¬ 
terest  for  two  reasons. 

First,  it  was  important  to  determine  if  a  full  account  of  the  circum¬ 
ferential  variation  of  the  wire  current,  due  to  the  edge  proximity,  would 
alter  the  input  impedance  as  determined  by  the  solution  in  section  B,  which 
did  not  account  for  any  variation.  As  will  be  «een,  the  solution  of  section 
B  yields  the  average  current  on  the  wire  and  hence  does  not  change  input 
impedance.  Secondly,  it  was  of  interest  to  determine  the  circumferential 
variation  of  the  current  at  the  base  of  the  wire  in  order  to  build  a  suit¬ 
able  expansion  mode  for  wires  attached  near  an  edge,  for  use  in  the  surface 
patch  MM  solution  of  Section  V. 

In  this  section  the  solution  of  section  B  is  extended  to  allow  for 
the  circumferential  variation  of  wire  current.  This  is  done  by  placing 
expansion  filaments  spaced  evenly  around  the  surface  of  the  wire,  i.e., 
by  constructing  a  wire-grid  model  of  the  wire  [91.  Each  filament  is  an 
independent  unknown,  and  can  also  be  broken  up  into  separate  modes  along 
its  length.  Test  modes  are  placed  alongside  the  expansion  filaments  spaced 
.OOOOOIX  away.  This  procedure  approximates  true  surface  testing  and  ex¬ 
pansion.  The  solution  of  section  B  is  then  a  special  case  (except  for 
test  mode  location)  of  this  more  general  formulation,  where  only  one  fila¬ 
ment  was  used. 

Figure  3a  shows  a  X/4  monopole,  radius  .OOlX,  placed  .05X  from  the 
edge  of  a  half-plane.  Also  shown  are  the  positions  of  eight  expansion 
filaments  used  in  this  example.  Below  is  Table  III,  showing  the  result¬ 
ing  filament  currents  at  the  base  of  the  monopole.  In  this  example  3  modes 
were  employed  in  the  longitudinal  direction  for  a  total  of  8x3  =  24  unknowns. 
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Figure  3.  Magnitude  and  phase  of  filament  current  at  the  base  of 
a  monopole  .05A  from  a  half-plane  edge,  vs.  position  around 
the  wire.  The  layout  of  the  eight  expansion  filaments 

is  also  shown. 


Table  III 


Filament  # 

Current  (ma) 

1 

3.1-j3.4 

2 

2.9-j2.5 

3 

2.3-j0.3 

4 

1.6+jl.9 

5 

1.2+j2.8 

6 

1.6+jl.9 

7 

2.3-j0.3 

8 

2.9-j2.5 

This  data  is  also  plotted  as  magnitude  and  phase  in  Figure  3b.  As 
expected,  the  data  show  symmetry  between  filaments  2,  3,  4  and  8,  7,  6, 
respectively.  The  data  also  show,  as  might  be  expected,  more  current  flow¬ 
ing  away  from  the  edge  than  toward  it.  If  the  filament  currents  in  the 
above  table  are  added  to  find  the  total  current,  the  "average"  input  im¬ 
pedance  can  be  calculated  and  is  found  to  be  identical  to  the  value  given 
by  the  simpler  one-filament  solution  of  section  8.  This  interesting  effect 
occurs  for  other  segmentation  schemes  and  distances  from  the  edge,  and 
is  similar  to  an  effect  reported  in  connection  with  closely  spaced  wire 
antennas  [91.  The  circumferential  variation  of  the  wire  current  becomes 
less  pronounced  as  one  moves  up  the  wire. 

E.  Expressions  for  Surface  Currents 

In  this  section  simple  closed  form  expressions  will  be  derived  for 
the  surface  current  density  on  a  perfectly  conducting  half-plane  induced 
by  an  infinitesimal  electric  dipole  on  the  surface  of  the  half-plane  and 
near  its  edge.  These  results  will  be  used  for  synthesis  of  an  expansion 
mode  for  a  wire  attachment  near  an  edge  in  the  surface  patch  MM  solution 
of  Section  V. 

The  surface  current  density  on  the  half-plane  is 

=  n  X  H  (16) 

For  the  total  surface  current  density,  one  must  vectorally  add  the  cur¬ 
rents  on  the  top  of  the  half-plane  (4>=0,  n=y)  to  those  on  the  bottom  of 
the  half-plane  (<t>=2iT,  n=-y): 

Jsx  =  Hz{4>=0)-H^(<))=2ti), 

(17) 

0s2  =  -H^(4)=0)+H^(^=27t). 
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lince  the  dipole  is  located  on  the  half-plane  (top)  at  2^=0,  ^^=0,  so  R=R 
and  Ir=Ir..  Then,  from  (3)  and  (4), 


TO  ^^0  I  I 

■^SX  ~  j  ■  3x  ^’^O  o 

I  °'<b=0  °  V2tt 


U  J  -2kz 


^  1 


Hp^kRj)  - 


°  ({.=0 


1  f^jPPol 

Now  let  a=sinh  "  - j,  then 

Mol  =  “  "  "PrI  ^ 

I  I  .  I  I  *  /Hp^kR  cosh  u)du  +  /“hp^kR  cosh  u)du 
*^*(t»=0  ^U=2ii  o  ^  ° 

-  T  HphkR  cosh  u)du  +  Hp^kR  cosh  u)du 
0  ^  0 

=  2  /* Hp^(kR  cosh  u)du 

'  0 

=  2^  (x,z|Xq,z^j)|^,^  ^0 


where 


4»(x,z|x-,z^)  = 


/  Hp^kR  cosh  u)du 
0 


Then  (18)  becomes 


-  ?k 


(21a) 


3 

3x 


i|»{x,z|x 


0*^0^ j 


J  -T^ 
sz  Tnj 


Hj^^kRj)  -  2k  1^  ,i,(x.z|Xq.Zq)' 


We  now  develop  a  small  argument  approximation  for  iji. 


(21b) 


For  small  x. 


H 


(2) 

1 


(X)  -  x/2  .  j  . 


SO  (20)  becomes 
<!»=  / 

0 


kR  cosh  u 
- 2 - 


+  .gJ 

iikR  cosh  u 


du 


1  ?  1  1  ^  J  1 

~  7  (4»"4>q)  :jj:^  tan  -  cos  j  (<J>-4)q) 


(22) 


For  the  present  case  4»=^q=0,  so  p=x  and  Pq=Xq.  As  a  check  on  the  accuracy 
of  this  approximation  ij»  was  calculated  by  numerical  integration  using  (20) 
and  compared  with  the  approximate  result  in  (22).  The  parameters  were 
x^=.05,  z=0,  .0075£X£.075,  X=1.,  and  the  data  are  listed  in  Table  IV  below: 


Table  IV 


X 

Exact  (Eq.  (20)) 

Approx.  (Eq.  (22)) 

.12+jl.90 

.12+jl.76 

.17+j3.20 

.17+j2.90 

.0225 

.21+j4.6I 

.21+j4.35 

.24+j6.98 

.24+j6.68 

.0375 

.27+jll.91 

.27+jll.57 

.29+j31.15 

.30+j30.76 

.0525 

.32+j63.08 

.32+j62.67 

. 34+ j 15. 39 

. 34+ j 14. 99 

.0675 

.36+j8.61 

.36+j8.23 

.0750 

.37+j5.92 

.38+j5.55 
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A  similar  check  was  made  after  was  differentiated  with  respect  to  x^. 

The  real  part  of  the  approximation  to  the  function  is  ignored  in  the 
remainder  of  the  analysis  since  it  is  negligible  in  comparison  to  the  imag- 
nary  part. 


The  last  step  is  to  carry  out  the  differentiations  as  indicated  in 

2 

(21).  When  this  is  done  there  will  be  present  terms  which  vary  as  1/R  . 
These  terms  are  due  to  the  charges  which  are  located  at  the  ends  of  the 
dipole.  Since  these  terms  vanish  in  a  physical  solution  which  has  con¬ 
tinuous  current,  they  are  omitted  in  the  final  expressions.  So,  for  a  ver 
tical  electric  dipole  of  moment  14  near  the  edge  of  a  half-plane,  the  sur¬ 
face  currents  are 


sz 


-I4z 


A/m^ 


(23) 


(24) 


Observe  that  with  Jx  dependence  and  l>(/x  dependence, 

as  x-K)  as  required  by  edge  conditions.  Also,  and  contain  the  cor¬ 
rect  1/R  singularity  in  the  vicinity  of  the  source.  Identifying  this  source 
singular  current  as  with  R  extending  radially  outward  from  the  source, 
gives 


^sR 


-214 

7^ 


XX 


-  A/m^, 


with  the  remainder  of  the  surface  current  as 


(25) 


It 

TT  Kj  J  0 

A/m^  (27) 

”  /rJ 
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contain  the 


Then  J  „  contains  the  source  singularity  and 

SK  SX 


edge 


and  J 


edge 


sz 


edge  terms.  Observe  that  is  dependent  on  circumferential  position 
around  the  source.  It  is  also  of  interest  to  note  that  Equations  (23)- 
(27)  are  not  frequency  dependent,  i.e.,  they  are  quasi-static  results. 
However,  retaining  the  first  term  in  the  real  part  or  higher  terms  in  the 
imaginary  part  of  (20)  will  introduce  a  frequency  dependence,  if  desired. 
Three-dimensional  plots  of  Equations  (23)-(24)  are  shown  in  Figures  4  and 
5. 


Note  that  the  surface  currents  are  unbounded  as  x^-K)  (as  the  source 
approaches  the  edge),  as  might  be  expected  on  the  basis  of  reciprocity. 
However,  if  the  magnitude  of  the  dipole  moment  is  decreased  at  the  same 
rate  as  the  currents  are  increasing  the  final  result  will  be  bounded  as 
x^-K).  This  allows  an  expression  for  surface  current  density  when  an  an¬ 
tenna  is  mounted  on  the  edge.  Thus,  replacing  the  dipole  moment  IH  by 
iJuXq  in  (23)  and  (24)  and  letting  x^+O  gives 


J 


SX 


(28) 


or. 


'sR 


-I 


/A 

N  X 


A/m 


(29) 


(30) 


with  _  Q_  Equation  (30)  is  the  form  of  the  surface  current 

density  for  a  vertical  dipole  mounted  right  at  the  edge  of  a  half-plane. 
Note  the  correct  source  and  edge  singularties.  A  3-D  plot  of  Equation 
(30)  is  shown  in  Figure  6. 
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III.  SOLUTION  FOR  MONOPOLE  NEAR  EDGE  OF  HALF-PLANE 
USING  MM/GTD 


A.  Formulation  of  MM/GTD  Solution 


In  this  section  the  solution  of-  section  II  for  a  vertical  monopole 
near  the  edge  of  a  half-plane  is  repeated  using  the  high  frequency  or  GTD 
Green's  function  for  the  diffracted  field  rather  than  the  exact  Green's 
function.  This  type  of  solution  is  termed  a  hybrid  MM/GTD  solution  [21. 


■’’he  geometrical  optics  (incident  and  reflected)  terms  are  separated 
from  the  total  field  and  used  to  compute  the  [Z^l  matrix,  as  in  section 
II.  The  fZl  matrix  is  computed  using  (14),  except  that  E^^^  is  now  cal¬ 
culated  using  the  Kouyoumjian-Pathak  GTD  formulation  for  spherical  wave 
incidence  [101  (using  the  notation  of  [10]): 

where,  for  our  geometry,  =  ti/2,  n=2. 


^  /V  _ 


IJ  =  z  z  D 


L  = 


PPr 


P+Pr 


-e 


-i  ^ 

J  4 


2j2^  sine'l 


F(kLa((|)-4,^))  F(kLa((j>+(i,Q))  | 


MX 

cosl 


,  2 


a(B)  =  2  cos'^  8/2,  F(x)  =  2jj7  /  e'^^  dj. 
Now,  e’  *2  =  0,  and  E^  =  E^  cos  <t),  so 
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(31) 


■E{(gE)  •  e-J"® 


In  (31)  eJ(Q^)  represents  the  incident  electric  field  evaluated  at  the 
half-plane  edge  which  is  transverse  to  the  ray  path. 


Because  PWS  modes  are  used  this  incident  field  can  be  expressed  in 
a  form  very  suitable  for  GTD  calculations.  In  particular,  the  total  exact 
electric  field  due  to  a  PWS  dipole  in  free-space  can  be  written  as  three 
spherical  waves  emanating  from  the  middle  and  endpoints  of  the  dipolefll]. 
These  spherical  waves  have  no  ray  path  field  components,  making  them  well 
suited  for  use  with  GTD.  Unfortunately,  this  decomposition  is  not  valid 
for  a  PWS  monopole,  i.e.,  a  ray  path  field  component  is  necessary  to  cor¬ 
rectly  express  the  field  from  a  PWS  monopole.  Thus,  the  incident  field 
from  the  bottom  mode  on  the  wire,  being  a  PWS  monopole,  is  not  evaluated 
exactly.  This  is  a  limitation  of  the  present  GTD. 


B.  Comparison  of  Results  from  Sections  II  and  III 

The  input  impedance  of  a  X/4  vertical  monopole  versus  distance  from 
the  edge  of  a  half-plane,  computed  using  MM/Green's  function  will  now  be 
compared  to  that  using  MM/GTD.  For  more  resolution,  the  change  in  input 
impedance,  AZ,  from  that  of  a  monopole  on  an  infinite  ground  plane  is  plotted. 
That  is,  AZ  =  Z^^-Z^,  where  Z^^  is  the  input  impedance  of  the  monopole 
on  the  half-plane  and  is  the  input  impedance  of  the  monopole  on  an  in¬ 
finite  ground  plane  (not  to  be  confused  with  the  [AZl  and  [Z^j]  matrices). 

The  impedance  Z^  and  matrix  [Z^l  is  the  same  regardless  of  which  method 
is  used. 


Figure  7  shows  AZ  vs.  d  for  a  X/4,  a=.004X  monopole  with  2  unknowns 
on  the  wire.  Note  that  the  MM/GTD  result  agrees  well  with  the  exact  MM/Green's 
function  result  for  the  monopole  as  close  as  .2X  from  the  edge.  The  good 
agreement  for  distances  greater  than  .2X  gives  confidence  in  the  results. 
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Figure  7a.  ReAZ  vs  d  for  MM/exact  Green's  function  and  for  MM/GTD 
solutions.  (^Z=Z.  -Z  ,  where  Z^  is  the  input  impedance  of  the 
vertical  V4  monopole  on  the  naif-plane  and  Z  is  the  input 
impedance  of  the  vertical  X/4  monopole  on  an 
infinite  ground  plane.) 
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Figure  7b.  Imff  vs  d  for  MM/exact  Green's  function  and  for 

MM/GTD  solutions. 
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As  mentioned  earlier,  the  MM/Green's  function  solution  becomes  un¬ 
stable  when  the  monopole  is  closer  than  about  .002X  from  the  edges.  The 
curves  in  Figure  7  were  thus  extrapolated  over  this  remaining  small  dis¬ 
tance,  but  compare  well  with  the  analysis  in  Section  IV  where  a  rigorous 
solution  is  done  for  a  monopole  mounted  on  the  edge  of  a  half-plane. 

IV.  SOLUTION  FOR  MONOPOLE  MOUNTED  AT  THE  EDGE  OF  A  HALF-PLANE 
USING  MM/GREEN'S  FUNCTION 

A.  Formulation  of  MM/Green's  Function  Solution 


The  geometry  of  the  problem  is  shown  in  Figure  8.  As  before,  the 
wire  and  half-plane  are  for  simplicity  taken  to  be  perfectly  conducting. 

The  monopole  lies  in  a  plane  perpendicular  to  the  edge  at  an  angle  a  from 
the  half-plane  surface,  connected  to  the  half-plane  at  the  coordinate  origin. 
Without  loss  of  generality,  a  is  restricted  to  0£o£ti. 

The  solution  is  similar  to  that  described  in  Section  II,  except  that 
here  the  exact  Green's  function  is  expressed  in  spherical  coordinates. 

The  test  modes  are  placed  on  the  wire  axis  (at  >*q.6q=  so 

the  required  T  field  component  to  compute  the  [Z1  matrix  is  E^^  (field 
in  r  direction  due  to  r  infinitesimal  dipole).  Note  that,  because  of  the 
non-zero  wire  radius,  the  field  along  the  wire  surface,  Ej^^,  actually  con¬ 
sists  of  both  an  E,.,.  and  an  Eq„  (field  in  6  direction  due  to  r  infinitesimal 
dipole)  component;  ^  -Eg^  ( -j,  where  p  is  the  cylindrical-co¬ 

ordinate  distance  of  the  field  point  from  the  edge  and  a  is  the  wire  radius. 
In  the  computation  of  fAZl  the  Eg^  term  is  ignored,  which  is  probably  justi- 
field  whenever  the  field  point  is  ten  wire  radii  or  more  from  the  edge 
(then  a/r«l).  It  is  felt  however,  that  ignoring  this  term  has  a  small 
but  noticeable  effect  on  the  input  impedance,  in  particular  on  the  imag¬ 
inary  part.  The  reason  that  this  Eg^  component  was  not  included  in  the 
analysis  was  that  it  was  too  difficult  to  do  so.  The  results  given  in 
the  next  section  show  that  this  approximation  is  justified. 
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HALF-PLANE 


Figure  8.  Geometry  and  coordinate  system  for  a  monopole  mounted 
at  the  edge  of  a  half-plane. 


From  another  paper  by  Senior [12],  the  radial  component  of  the  electric 
field  due  to  a  radially  directed  electric  dipole,  with  moment  14,  near 
the  edge  of  a  half-plane  is 


0  J  -J 


The  magnetic  field  components  tangential  to  the  half-plane,  for  use  later 
in  determining  the  surface  current,  are 


n 


0r  4Ttj 


=  0 
rr 


_T  )/■ 

r  sine  [  ^R'M 


(33) 


where  and  are  defined  in  (5).  Equations  (32),  (33)  represent  the 
fields  in  terms  of  a  Debye  potential,  whereas  Equations  (2)-(4)  represented 
the  fields  in  terms  of  a  Hertz  potential. 


Now,  it  can  be  shown  that 


e-jkR 

m — 


1 

i 


and 


r 


(34) 


(35) 


where  the  minus  sign  in  (35)  is  due  to  the  fact  that  the  image  dipole  is 
in  the  -r  direction.  Then,  defining 


c  _  rinc  ^  .image  ^  .dif 
rr  rr  rr  rr 

we  have 


pdif 

*^rr  4irj 

0 

(36) 


(37) 


dp  dR  * 

where  tt  and  ti  are  defined  in  (10)  or  (11).  Observe  that  when  a-"T/2 
there  is  no  reflected  E^^  field  for  4»Tr/2,  so  (37)  does  not  strictly  repre¬ 
sent  the  diffracted  field  in  the  usual  sense. 

The  source  and  image  fields  are  then  separated  from  the  total  field, 
and  the  [Z^l  section  of  the  total  impedance  matrix  [Z]  is  computed  in  closed 
form  as  before,  using  Richmond's  subroutines  [81.  The  (AZl  elements  are 
found  by 
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'  ■  /  /  (r|r^)Ij(r)dr  dr„,  (38) 

*"01  '"j 

where  r^.  ranges  over  the  test  mode  and  rj  ranges  over  the  expansion  mode. 

The  differentiations  in  (37)  can  be  eliminated  by  integration  by  parts. 

r 

If  we  let  KCrlr^)  =  (Ip-Ij^,),  much  simplification  is  possible  by  noting 

that  K(r|0)=0,  K(r|rQ)|^  .=0,  and  ^  I(r)=-k2l(r)-2kcot  kdafr-r*") 

for  I(r)  being  the  current  of  a  PWS  dipole  and  r'"  the  midpoint  (terminal 
point)  of  the  dipole.  The  final  form  after  simplification  is 


_ 

sin  kd 


4iij 


/I  .(r)  [K(r|rJ.)+K(r|rQ.)-2cos  kd  K(r|rQ-)]dr, 

(39) 


where  r^^,  r^^,  r^.  are  the  radial  coordinates  of  the  top,  middle,  bottom, 
respectively  of  the  i-th  test  dipole.  If  the  test  mode  is  the  bottom  mode, 

_n  ^ 

which  is  a  monopole,  then  r^.  =  r°.  =0.  d  is  the  length  of  the  mode  seg¬ 
ments. 


This  solution  was  programmed  and  run  for  aX/4,  a=.00lA  monopole, 
with  a  delta  gap  generator  at  the  base  of  the  monopole,  vs  angle  a. 


Below  are  shown  the  AZ^.j  elements  with  two  wire  modes  and  a=ir/2. 

The  fZ^l  matrix  is  identical  to  that  given  in  Section  II.  The  expansion 
modes  are  numbered  from  top  to  bottom  (attachment  point)  as  befor  ;. 


[AZ] 


(.75-j9.56)  (7.45+j5.03) 

(7.43-j24.29)  (13.38+j23.88) 


The  asymmetry  in  the  [AZ]  matrix  is  attributed  to  the  fact  that  the 

“r 

field  component  was  neglected,  and  possibly  because  full  surface  testing 
and  expansion  was  not  used.  However,  the  lAZ]  elements  are  relatively 
small  compared  to  the  IZq]  elements  (see  section  IIC)  reducing  the  final 
error. 
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Figure  13  shows  using  two  modes  on  the  wire  compared  with  the 
surface  patch  MM  solution.  This  result  at  a=TT/2  also  agrees  with  the  analysis 
for  the  vertical  monopole  in  Section  II  as  approaches  zero.  Note  that 

^in  goes  to  zero  as  a  approaches  zero  as  expected,  since  the  antenna  then 
becomes  shorted  by  the  ground  plane.  As  a  approaches  zero,  the  reactance 
curve  takes  a  sharp  turn  downward  before  going  to  zero.  This  effect  is 
probably  explained  by  the  fact  that  the  capacitance  between  the  monopole 
and  half-plane  increases  sharply  as  the  monopole  gets  closer  to  the  half¬ 
plane,  before  going  to  zero  as  it  touches  the  half-plane. 

Another  independent  check  is  available  from  Tai[5]  who  computed  the 
input  resistance  of  a  monopole  at  a=n  by  integrating  the  far  field  radi¬ 
ation  pattern.  Tai  obtained,  using  one  mode  on  the  wire,  a  value  of  86.3;2 
while  our  analysis,  using  one  mode,  gives  85.2il  for  the  input  resistance 
(the  input  reactance  cannot  be  calculated  using  Tai's  method).  Also  shown 
in  Figure  13  is  a  surface  patch  result  obtained  by  Glisson(13]  for  a  X/4 
monopole  mounted  on  the  edge  of  a  finite  plate,  which  compares  reasonably 
with  the  other  results. 

B.  Expression  for  Surface  Current 

In  this  section  a  simple  closed  form  expression  is  derived  for  the 
surface  current  density  on  a  perfectly  conducting  half-plane  induced  by 
an  infinitesimal  electric  dipole  mounted  at  the  edge.  The  surface  current 
density  on  the  half-plane  is  =  n  x  H.  Since  n  =  -4)  and  =  0  the  only 
surface  current  component  is  =  -Hg^((j)=0)+Hg^(4i=2ir) ,  which  is  the  sum 
of  the  surface  currents  on  the  top  and  bottom  of  the  half-plane.  As  ex¬ 
pected,  the  entire  surface  current  flows  radially  outward  from  the  dipole. 
Then,  using  (33), 

-we  ^^^0  3 

"^sr  “  2r^sTn0  Aitj  3^  ^  l(j>=0  ^R"^R' ^ 
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Now,  R(<t)=0)=R((t»=2TT)  and  ,  so  the  /  (  )  part  of  the  in- 

tegrals  in  Ip  and  Ip,  of  (5)  cancel,  since  they  are  only  functions  of  R 
or  R*.  Also,  it  can  be  shown  that  the  remaining  terms  can  be  written  as 

1  -  H  I 

sr  r^sTHS  I  ,^0  ’ 

where  t  is  defined  in  (20).  Then,  using  the  small  argument  approximation 
given  in  (22), 


_  -4ja)e 


sr  ~  Tfkr^sine  4iTj 


■0  3  Z'  1  *  -1  [^J^o  1 


cos  ^  (4.-4)^ 


Carrying  out  the  differentiation  and  ignoring  terms  of  order  1/R^, 
which  account  for  end  charges,  gives 

-4J»e  ‘“o  1  f  *0  .  .  .  .  r-  .  *0 

’  " — ^ -  —  — '  2 -  cos sin<p^smd+Jp^  sin 


sr  “  itkr^sine  4Tij  p2 


Up  to  this  point  the  infinitesimal  dipole  has  been  located  at  ('"q/0, 
0Q=Tt/2,4)Q) .  To  move  the  dipole  to  the  edge  (>"^-*0)  the  dipole  moment  must 
be  reduced  as  Jr^,  as  was  done  in  Section  II,  to  maintain  finite  fields. 
Thus,  for  a  dipole  of  moment  l/I»^  at  the  edge,  the  surface  current  density 
is  (Rj=r  at  r^=0) 

-ijlo  1  *^o 

J-_  =  -5 — —  -y  sin  -5-  ,  or,  since  sinG  =  p/r  =  x/r, 
n  sine  r'^  ^ 


•^sr  =  7  yi  *  ?  5’" 


p-  A/m 


For  a  vertical  dipole,  giving 
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Equation  (45)  agrees  in  form  with  result  (30)  in  Section  II  but  differs  f 

by  a  factor  of  J?.  The  reason  for  this  discrepancy  is  not  known  at  present,  f 

f 

but  probably  has  something  to  do  with  the  strength  of  the  dipole  moment. 

I 

For  this  work,  we  are  primarily  interested  in  the  form  of  the  current.  j; 

Observe  that  Equations  (44)-(45)  have  the  proper  source  and  edge  singu-  f 

larities,  and  that  the  orientation  angle  of  the  dipole,  (J)^,  does  not  affect 

0  ,  t 

the  form  of  the  current  density. 

Glisson  (13]  derived  a  similar  result  for  the  surface  currents  in¬ 
duced  by  a  radial  source  mounted  on  the  edge  of  a  wedge  of  arbitrary  angle.  ; 

For  the  case  of  a  half-plane  his  result  reduces  to  J  _  ,  which  disagrees 

/xr  If 

with  the  form  of  our  result  in  terms  of  the  source  singularity.  It  is  j: 

felt  that  Glisson's  solution  does  not  predict  the  correct  source  behavior 
because  his  solution  was  based  on  retaining  only  the  first  term  in  an  eigen-  i! 

function  expansion.  As  mentioned  previously,  this  procedure  usually  pre-  ' 

diets  the  correct  edge  behavior  but  not  the  correct  source  behavior.  | 

1 

I 

V.  SURFACE  PATCH  MM  SOLUTION 

This  section  deals  with  methods  of  modifying  the  author's  surface  patch 
modelling  program  (41  to  handle  antennas  mounted  near  or  at  the  edge  of 
plates.  Originally,  the  program  used  an  attachment  mode  which  had  the 
following  features:  \ 

i 

a)  enforced  continuity  of  current  from  wire  to  plate.  : 

I 

I 

b)  enforced  the  proper  p/p  dependence  for  surface  current  in  the  i 

vincinity  of  the  attachment  point. 


d)  worked  well  for  wire  attachments  .IX  or  further  from  an  edge. 


This  attachment  mode,  shown  in  Figure  9a,  failed  for  antennas  mounted 
closer  than  .lA  from  the  edge  as  shown  in  Figures  11  and  12,  probably  for 
the  following  reasons: 

a)  disk  mode  radius  was  limited  to  the  distance  of  the  monopole 
from  the  edge,  thus  reducing  the  area  in  which  the  p/p  singu¬ 
larity  was  enforced. 

b)  disk  mode  did  not  contain  any  (j)-variation  in  current  density, 
which  probably  becomes  increasingly  important  as  the  antenna 
gets  closer  to  the  edge. 

c)  for  small  distances  from  the  edge,  the  current  filaments  con¬ 
stituting  the  disk  mode  are  very  short,  making  calculation  of 
mutual  impedances  difficult. 

Thus,  a  better  attachment  mode  was  needed  to  handle  antennas  mounted  close  to 
or  at  the  edge  of  plates. 

A.  Description  of  Attachment  Mode  Development 

All  of  the  attachment  modes  that  were  developed  consisted  of  a  disk 
section  which  rested  on  the  plate  structure,  and  a  vertical  wire  monopole 
with  terminals  at  the  center  of  the  disk.  Since  this  wire  monopole  was 
common  to  all  attachment  modes,  it  is  not  described  in  detail  here. 

Mode  A)  The  original  disk  monopole  part  of  the  attachment  mode  had  surface 
current  density  [4] 


where  a  is  the  wire  radius  and  b  is  the  outer  radius  of  the  disk,  usually 
taken  to  be  between  O.IX  and  0.25X,  if  room  permits  (see  Figure  3  of  [4]). 
Figure  9a  shows  how  this  mode  would  be  applied  to  a  wire  attachment  near 
an  edge. 


Mode  B)  The  first  approach  toward  improvement  of  the  disk  mode  A)  was 
to  maintain  the  disk  radius  at  about  .2X  where  it  overlapped  the  plate, 
but  cut  off  the  disk  where  it  overhung  the  plate,  see  Figure  9b.  The  surf- 
face  current  density  for  this  disk  mode  can  be  written  as 


Jp(p.<J»)  =  B 


sin  k(b(4))-p) 

2  up 


(47) 


where  B  is  a  constant  chosen  for  unit  terminal  current,  i.e., 

2ti 

/  J  (a,(D)ad4.  =  lA,  (48) 

(i)=0  P 

and  b((t))  is  the  outer  radius  of  the  disk  mode,  now  expressed  as  a  non-con¬ 
stant  function  of  (|). 


This  type  of  mode  overcomes  some  of  the  drawbacks  of  the  previous 
mode  in  that  it  has  a  (j)-vari ation  of  current  and  yields  a  larger  area  where 
the  1/p  singularity  is  enforced.  The  results  obtained  when  using  this 
mode  are  not  significantly  better  than  those  obtained  with  the  original 
mode  A).  This  is  probably  due  to  the  fact  that  the  ({(-dependence  is  not 
close  enough  to  the  exact  dependence.  For  instance,  when  p  is  close  to 
a,  the  current  density  of  mode  B)  is  large  but  does  not  vary  appreciably 
with  ({),  while  the  analysis  of  Section  II,  Equation  (25),  shows  that  sig¬ 
nificantly  less  current  is  flowing  toward  the  edge  than  away  from  it.  An¬ 
other  point  is  that  this  mode  does  not  have  a  uniform  limit  as  the  v’ttach- 
ment  point  approaches  the  edge.  One  would  expect  that  the  current  flowing 
toward  the  edge  would  smoothly  decrease  to  zero  as  the  attachment  moves 
to  the  edge.  Mode  B)  does  not  have  this  behavior.  The  disk  modes  pre¬ 
sented  below  do  have  this  desired  asymmetry. 
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Mode  C)  The  next  step  was  to  construct  a  disk  mode  which  had  less  current 
flowing  toward  the  edge  than  away  from  it.  One  distribution  which 
does  this  is 

=  C  — sin  kb((|)),  (49) 

where  C  is  a  constant  chosen  to  satisfy  (48).  This  current  distribution 
roughly  corresponds  to  that  obtained  from  the  exact  solution  for  the  cir¬ 
cumferential  current  variation,  as  shown  in  Figure  3b.  This  type  of  mode 
has  the  desired  property  of  uniformly  approaching  a  semi-circular  disk 
as  the  attachment  point  moves  to  the  edge  (Figure  9c).  Also,  this  mode 
reduces  to  the  standard  attachment  mode,  mode  A),  when  the  distance  be¬ 
tween  the  attachment  point  and  edge  becomes  greater  than  the  disk  outer 
radius  (_.2X),  since  b(4i)  =  constant  for  this  case. 

This  type  of  attachment  mode  gave  the  best  results  for  antennas  at 
or  near  an  edge.  These  results  are  compared  with  others  later. 

Mode  D)  The  next  type  of  disk  mode  which  was  developed  was  one  constructed 
from  an  arbitrary  number  of  wedge-shaped  sectors  of  current,  with  current 
distribution  like  that  of  disk  mode  B)  over  each  sector,  but  with  each 
sector  being  an  independent  unknown.  See  Figure  10,  which  shows  a  4  sec¬ 
tor  disk.  This  type  of  mode  has  the  advantage  of  allowing  the  MM  solution 
to  arrive  at  currents  (both  magnitude  and  phase)  which  best  satisfy  the 
integral  equation  in  a  least  mean  square  sense.  The  resolution  of  the 
variation  in  disk  current  can  be  made  finer  by  increasing  the  number  of 
wedge-shaped  sectors  that  the  disk  is  broken  into.  When  eight  such  sec¬ 
tors  were  used  the  results  for  input  impedance  agreed  fairly  well  with 
those  obtained  using  mode  C).  The  resulting  sector  currents  for  x^=.05X, 
a=.001X,  and  8  sectors  with  sector  1  oriented  directly  away  from  the  edge  are 
given  in  Table  V  below: 
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Figure  9. 


Note  the  symmetry  between  sectors  (4,6),  (3,7)  and  (2,8);  the  sym¬ 
metry  is  not  exact  probably  because  of  numerical  integration  error.  This 
current  distribution  is  similar  in  magnitude  to  that  shown  in  Figure  3b, 
but  not  in  phase.  This  is  interesting,  because  it  implies  that  the  en¬ 
forcement  of  the  proper  phase  variation  is  not  crucial  to  a  correct  result 
for  input  impedance.  This  could  explain  why  mode  C),  which  enforces  cur¬ 
rent  magnitude  but  not  phase,  gives  good  results. 

Edge  Mode: 


Recently,  there  has  been  an  interest  in  directly  enforcing  the  edge 
condition  in  MM  solutions  [14].  Thus,  it  was  decided  to  construct  a  mode, 
based  on  the  analysis  of  Section  HE,  that  would  enforce  the  edge  singu¬ 
larity  of  the  surface  current  component  parallel  to  the  edge  of  the  half¬ 
plane,  and  study  its  effect  on  the  solution.  This  "edge  mode"  was  con¬ 
structed  from  PWS  current  filaments  to  embody  the  essential  features  of 
Equation  (27).  The  surface  current  density  of  this  mode  is 


J5(x,z) 


A 


r 

-sin  k(z+g.) 
sin  k(?--ZQ) 

sin  kz 

1  sin  kz- 


sin  k(z-jl) 

sin  k(z'-£y 


for  -8,<z<-Zq, 
for  -Zq<z<Zjj 

for  Zq<z<JI 


(50) 


where  the  coordinate  system  is  shown  in  Figure  1,  A  is  a  suitable  normal¬ 
izing  constant,  and  2a  is  the  overall  length  of  the  mode.  The  width  of 
the  mode  (extent  in  x-direction)  is  arbitrary,  but  since  most  of  the  cur¬ 
rent  is  located  near  the  singularity  at  x=0,  this  width  is  not  too  critical 
if  chosen  to  be  .05X  or  greater.  The  result  of  using  this  edge  mode  with 
different  lengths  and  widths  was  interesting.  First,  the  edge  mode  was 
excited  in  the  surface  patch  MM  solution,  as  indicated  by  a  non-zero  mode 
current.  However,  the  input  impedance  did  not  significantly  change  from 
the  result  obtained  without  using  an  edge  mode.  Thus,  all  of  the  other 
surface  patch,  wire,  and  attachment  currents  changed  in  the  presence  of 
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the  edge  mode  in  such  a  way  as  to  leave  the  input  impedance  unaltered. 

At  the  present  time  the  best  explanation  we  have  for  this  effect  is  that 
the  singular  behavior  of  the  edge  current  can  be  approximated  adequately 
in  a  least-mean-square  sense  by  the  surface-patch  dipoles  alone.  Incor¬ 
poration  of  a  special  edge  mode  explicitly  satisfying  the  edge  condition 
does  not  improve  on  this  situation  in  the  present  case. 

B.  Results  and  Comparisons  with  Solutions 
from  Sections  II.  Ill,  IV 

In  this  section  results  from  other  sections  are  compared  with  the 
surface  patch  MM  results.  Of  necessity,  the  surface  patch  MM  solution 
must  be  based  on  a  finite  size  plate  rather  than  a  semi-infinite  half-plane. 
For  this  reason  it  is  more  appropriate  to  plot  rather  than  AZ.  Also, 
the  size  of  this  finite  ground  plane  can  have  some  effect  on  the  input 
impedance,  as  shown  in  Table  VI  below  which  gives  the  input  impedance  of 
a  X/4  monopole  .05X  away  from  the  long  edge  of  different  sized  ground  planes. 

Table  VI 


Plate  Size 

^in 

.4 

X.5A 

44.+jl3. 

.8 

xl.X 

71.+jl2. 

.8 

xl.2X 

54.-j4.6 

.8xl.5A 

42+jl2. 

This  variation  seems  to  be  a  plate  resonance  effect,  and  probably 
explains  the  offset  of  a  few  ohms  between  the  surface  patch  MM  and  other 
solutions  in  the  following  figures. 

Figures  11  and  12  show  the  real  and  imaginary  parts  of  the  input  im¬ 
pedance  of  a  A/4,  a=.00lA  monopole  at  a  distance  d  from  the  edge  of  a  half 
plane,  as  calculated  with  the  exact  MM/Green's  function  solution  (Section 


jure  11.  Input  resistance  of  a  X/4  monopole  vs  distance,  d,  from  the  edge  of  a  half-plane 
computed  using  MM/exact  Green's  function,  MM/6TD,  surface  patch  with  old  att.  mode  (A), 

and  surface  patch  with  new  att.  mode  (C). 


INPUT  IMPEDANCE  OF  MONOPOLE  (  X /4 ) 
NEAR  EDGE  OF  HALF  -  PLANE 


L. 


II),  the  MM/GTD  solution  (Section  III),  and  the  surface  patch  MM  solution 
using  mode  A)  and  mode  C).  Remember  that  the  surface  patch  MM  results 
are  for  finite  ground  planes,  in  this  case  .4Xx.5X,  with  7  plate  modes, 

1  attachment  mode,  and  1  wire  mode.  The  solutions  of  Sections  II  and  III 
both  used  2  modes  on  the  wire  monopole.  From  these  graphs  it  can  be  seen 
that  the  original  surface  patch  attachment  mode  does  indeed  fail  when  the 
antenna  is  chosen  closer  than  .IX  to  the  edge,  and  that  the  modified  at¬ 
tachment  mode  (mode  C)  performs  much  better.  Recall  that  the  offset  of 
the  surface  patch  results  from  the  other  curves  can  be  accounted  for  by 
the  finite  plate  size  required  for  the  surface  patch  result. 

Figure  13  presents  the  results  for  a  X/4,  a=.OOlX  monopole  mounted 
at  the  edge  of  a  half-plane  vs  angle  based  on  the  solution  discussed  in 
Section  IV,  and  compared  with  the  surface  patch  MM  solution  with  attachment 
mode  C).  Discussion  of  the  solution  of  Section  IV  and  comparison  with 
independent  results  was  given  in  Section  IV.  The  surface  patch  MM  result 
in  Figure  13  used  a  semicircular  attachment  mode  like  that  shown  in  Figure 
9c,  with  a  .4xx.5x  plate  and  9  modes  total. 

VI.  CONCLUSION 

This  report  has  presented  analysis  and  results  of  the  problem  of  a 
monopole  antenna  mounted  near  or  at  the  edge  of  a  half -plane  or  plate. 

The  canonical  solutions  using  exact  Green's  functions  have  been  presented 
as  a  basis  for  comparison  with  other  solutions.  Expressions  for  near-zone 
surface  current  density  were  derived  and  used  to  synthesize  an  attachment 
mode  to  handle  wire  attachments  near  plate  edges  in  a  general  surface  patch 
MM  solution. 

The  most  obvious  extension  of  the  work  presented  here  would  be  the 
analysis  of  monopoles  mounted  near  or  on  the  edge  of  a  wedge  of  arbitrary 
angle.  The  handling  of  wedges  in  the  canonical  solutions  would  be  more 
difficult  than  that  for  half-planes  because  the  Green's  function  for  a 
wedge  is  not  in  as  convenient  a  form  as  is  the  half-plane  Green's  function. 
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INPUT  IMPEDANCE  OF  X/4  MONOPOLE 
AT  EDGE  OF  HALF- PLANE  vs.  ANGLE 


Figure  13.  Input  impedance  Qf  a  X/4  radial  monopole  mounted  at  the  edge 
of  a  half-plane  vs.  angle,  computed  using  MM/exact  Green's  function 
and  surface  patch  with  att.  mode  (C).  Also  shown  are  independent 
results  for  R  and  Z  at  a=180°,  from  Tai  and  Glisson,  respectively. 


I 
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However,  the  surface  patch  solution  for  antennas  would  probably  be  easier, 
for  two  reasons.  First,  the  half-plane  represents  the  "worst  case"  as 
far  as  edge  effect.  That  is,  a  wedge  with  interior  angle  greater  than 
zero  would  be  expected  to  have  less  effect  on  input  impedance  than  would 
a  half-plane.  Second,  the  surface  current  flow  around  the  edge  of  a  wedge 
can  probably  be  handled  with  surface  patch  modes  and  minor  modifications 
to  the  attachment  mode. 

Another  related  area  for  future  work  would  be  the  case  of  a  monopole 
mounted  near  a  corner  or  vertex.  Also,  these  solutions  could  be  extended 
to  account  for  mutual  coupling  between  two  or  more  antennas  near  edges 
in  a  straightforward  manner. 
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